Mismatch Repair Protein Msh2 Is Necessary for Macronuclear Stability and Micronuclear Division in Tetrahymena thermophila.

Mismatch repair (MMR) is a conserved mechanism that is primarily responsible for the repair of DNA mismatches during DNA replication. Msh2 forms MutS heterodimer complexes that initiate the MMR in eukaryotes. The function of Msh2 is less clear under different chromatin structures. Tetrahymena thermophila contains a transcriptionally active macronucleus (MAC) and a transcriptionally silent micronucleus (MIC) in the same cytoplasm. Msh2 is localized in the MAC and MIC during vegetative growth. Msh2 is localized in the perinuclear region around the MIC and forms a spindle-like structure as the MIC divides. During the early conjugation stage, Msh2 is localized in the MIC and disappears from the parental MAC. Msh2 is localized in the new MAC and new MIC during the late conjugation stage. Msh2 also forms a spindle-like structure with a meiotic MIC and mitotic gametic nucleus. MSH2 knockdown inhibits the division of MAC and MIC during vegetative growth and affects cellular proliferation. MSH2 knockdown mutants are sensitive to cisplatin treatment. MSH2 knockdown also affects micronuclear meiosis and gametogenesis during sexual development. Furthermore, Msh2 interacts with MMR-dependent and MMR-independent factors. Therefore, Msh2 is necessary for macronuclear stability, as well as micronuclear mitosis and meiosis in Tetrahymena.

Amitosis as a strategy of cell division-Insight from the proliferation of Tetrahymena thermophila macronuclei.

Cell division is a necessity of life which can be either mitotic or amitotic. While both are fundamental, amitosis is sometimes considered a relic of little importance in biology. Nevertheless, eukaryotes often have polyploid cells, including cancer cells, which may divide amitotically. To understand how amitosis ensures the completion of cell division, we turn to the macronuclei of ciliates. The grand scheme governing the proliferation of the macronuclei of ciliate cells, which involves chromosomal replication and amitosis, is currently unknown, which is crucial for developing population genetics model of ciliate populations. Using a novel model that encompasses a wide range of mechanisms together with experimental data of the composition of mating types at different stages derived from a single karyonide of Tetrahymena thermophila, we show that the chromosomal replication of the macronucleus has a strong head-start effect, with only about five copies of chromosomes replicated at a time and persistent reuse of the chromosomes involved in the early replication. Furthermore the fission of a fully grown macronucleus is non-random with regard to chromosome composition, with a strong tendency to push chromosomes and their replications to the same daughter cell.

An Extremely Streamlined Macronuclear Genome in the Free-Living Protozoan Fabrea salina.

Ciliated protists are among the oldest unicellular organisms with a heterotrophic lifestyle and share a common ancestor with Plantae. Unlike any other eukaryotes, there are two distinct nuclei in ciliates with separate germline and somatic cell functions. Here, we assembled a near-complete macronuclear genome of Fabrea salina, which belongs to one of the oldest clades of ciliates. Its extremely minimized genome (18.35 Mb) is the smallest among all free-living heterotrophic eukaryotes and exhibits typical streamlined genomic features, including high gene density, tiny introns, and shrinkage of gene paralogs. Gene families involved in hypersaline stress resistance, DNA replication proteins, and mitochondrial biogenesis are expanded, and the accumulation of phosphatidic acid may play an important role in resistance to high osmotic pressure. We further investigated the morphological and transcriptomic changes in the macronucleus during sexual reproduction and highlighted the potential contribution of macronuclear residuals to this process. We believe that the minimized genome generated in this study provides novel insights into the genome streamlining theory and will be an ideal model to study the evolution of eukaryotic heterotrophs.

Micromanipulation in Paramecium: From non-mendelian inheritance to the outlook for versatile micromachines.

This review addresses nine areas of knowledge revealed by micromanipulations performed with Paramecium. Microinjection has shown that sexual maturation and senescence of Paramecium caudatum is a programmed process conducted by a specific gene and its product protein. In Paramecium tetraurelia, autogamy was revealed to depend on the number of DNA syntheses rather than the number of cell divisions in clonal aging. The cytoplasmic complementarity test established that microinjection of wild-type cytoplasm can correct genetic defects of mutants. The concept of complementarity together with protein chemistry revealed compounds that control membrane excitability. In non-Mendelian inheritance, noncoding small RNAs made from the parental micronucleus regulate the rearrangement of the progeny's macronuclear DNA. The macronucleus has the potential to be used as a factory for genetic engineering. The development and differentiation of progeny's nuclei in mating pairs are controlled by the parental macronucleus. The chemical reaction processes associated with exocytosis have been revealed by microinjection of various enzymes and antibodies. Using the fusion gene of histone H2B and yellow-fluorescence protein, it was revealed that the fusion gene-mRNA is transferred between cells during mating. Experiments with endosymbiotic bacteria and the host shed light on the conditions needed to establish sustainable symbiotic relationships.

Absolute quantification of chromosome copy numbers in the polyploid macronucleus of Tetrahymena thermophila at the single-cell level.

Amitosis is widespread among eukaryotes, but the underlying mechanisms are poorly understood. The polyploid macronucleus (MAC) of unicellular ciliates divides by amitosis, making ciliates a potentially valuable model system to study this process. However, a method to accurately quantify the copy number of MAC chromosomes has not yet been established. Here, we used droplet digital PCR (ddPCR) to quantify the absolute copy number of the MAC chromosomes in Tetrahymena thermophila. We first confirmed that ddPCR is a sensitive and reproducible method to determine accurate chromosome copy numbers at the single-cell level. We then used ddPCR to determine the copy number of different MAC chromosomes by analyzing individual T. thermophila cells in the G1 and the amitotic (AM) phases. The average copy number of MAC chromosomes was 90.9 at G1 phase, approximately half the number at AM phase (189.8). The copy number of each MAC chromosome varied among individual cells in G1 phase and correlated with cell size, suggesting that amitosis accompanied by unequal cytokinesis causes copy number variability. Furthermore, the fact that MAC chromosome copy number is less variable among AM-phase cells suggests that the copy number is standardized by regulating DNA replication. We also demonstrated that copy numbers differ among different MAC chromosomes and that interchromosomal variations in copy number are consistent across individual cells. Our findings demonstrate that ddPCR can be used to model amitosis in T. thermophila and possibly in other ciliates.

Early developmental, meiosis-specific proteins - Spo11, Msh4-1, and Msh5 - Affect subsequent genome reorganization in Paramecium tetraurelia.

Developmental DNA elimination in Paramecium tetraurelia occurs through a trans-nuclear comparison of the genomes of two distinct types of nuclei: the germline micronucleus (MIC) and the somatic macronucleus (MAC). During sexual reproduction, which starts with meiosis of the germline nuclei, MIC-limited sequences including Internal Eliminated Sequences (IESs) and transposons are eliminated from the developing MAC in a process guided by noncoding RNAs (scnRNAs and iesRNAs). However, our current understanding of this mechanism is still very limited. Therefore, studying both genetic and epigenetic aspects of these processes is a crucial step to understand this phenomenon in more detail. Here, we describe the involvement of homologs of classical meiotic proteins, Spo11, Msh4-1, and Msh5 in this phenomenon. Based on our analyses, we propose that proper functioning of Spo11, Msh4-1, and Msh5 during Paramecium sexual reproduction are necessary for genome reorganization and viable progeny. Also, we show that double-strand breaks (DSBs) in DNA induced during meiosis by Spo11 are crucial for proper IESs excision. In summary, our investigations show that early sexual reproduction processes may significantly influence later somatic genome integrity.

Macronuclear development in ciliates, with a focus on nuclear architecture.

Ciliates are defined by the presence of dimorphic nuclei as they have both a somatic macronucleus and germline micronucleus within each individual cell. The size and structure of both germline micronuclei and somatic macronuclei vary tremendously among ciliates. Except just after conjugation (i.e. the nuclear exchange in their life cycle), the germline micronucleus is transcriptionally inactive and contains canonical chromosomes that will be inherited between generations. In contrast, the transcriptionally active macronucleus contains chromosomes that vary in size in different classes of ciliates, with some lineages having extensively fragmented gene-sized somatic chromosomes while others contain longer multigene chromosomes. Here, we describe the variation in somatic macronuclear architecture in lineages sampled across the ciliate tree of life, specifically focusing on lineages with extensively fragmented chromosomes (e.g. the classes Phyllopharyngea and Spirotrichea). Further, we synthesize information from the literature on the development of ciliate macronuclei, focusing on changes in nuclear architecture throughout life cycles. These data highlight the tremendous diversity among ciliate nuclear cycles, extend our understanding of patterns of genome evolution, and provide insight into different germline and somatic nuclear features (e.g. nuclear structure and development) among eukaryotes.

The hidden genomic diversity of ciliated protists revealed by single-cell genome sequencing.

BACKGROUND: Ciliated protists are a widely distributed, morphologically diverse, and genetically heterogeneous group of unicellular organisms, usually known for containing two types of nuclei: a transcribed polyploid macronucleus involved in gene expression and a silent diploid micronucleus responsible for transmission of genetic material during sexual reproduction and generation of the macronucleus. Although studies in a few species of culturable ciliated protists have revealed the highly dynamic nature of replicative and recombination events relating the micronucleus to the macronucleus, the broader understanding of the genomic diversity of ciliated protists, as well as their phylogenetic relationships and metabolic potential, has been hampered by the inability to culture numerous other species under laboratory conditions, as well as the presence of symbiotic bacteria and microalgae which provide a challenge for current sequencing technologies. Here, we optimized single-cell sequencing methods and associated data analyses, to effectively remove contamination by commensal bacteria, and generated high-quality genomes for a number of Euplotia species. RESULTS: We obtained eight high-quality Euplotia genomes by using single-cell genome sequencing techniques. The genomes have high genomic completeness, with sizes between 68 and 125 M and gene numbers between 14K and 25K. Through comparative genomic analysis, we found that there are a large number of gene expansion events in Euplotia genomes, and these expansions are closely related to the phenotypic evolution and specific environmental adaptations of individual species. We further found four distinct subgroups in the genus Euplotes, which exhibited considerable genetic distance and relative lack of conserved genomic syntenies. Comparative genomic analyses of Uronychia and its relatives revealed significant gene expansion associated with the ciliary movement machinery, which may be related to the unique and strong swimming ability. CONCLUSIONS: We employed single-cell genomics to obtain eight ciliate genomes, characterized the underestimated genomic diversity of Euplotia, and determined the divergence time of representative species in this subclass for the first time. We also further investigated the extensive duplication events associated with speciation and environmental adaptation. This study provides a unique and valuable resource for understanding the evolutionary history and genetic diversity of ciliates.

[Packaging of Subchromosomal-Size DNA Molecules in Chromatin Bodies in the Ciliate Macronucleus].

A fundamental difference between somatic nuclei (macronuclei) of ciliates and cell nuclei of higher eukaryotes is that the macronuclear genome is a huge number (up to tens or hundreds of thousands) of gene-sized (0.5-25 kb) or subchromosomal (up to 2000 kb) minichromosomes. Electron microscopy shows that macronuclear chromatin usually looks like chromatin bodies or fibrils 200-300 nm thick in the interphase. However, the question of how many DNA molecules are contained in an individual chromatin body remains open. The organization of chromatin in macronuclei was studied in the ciliates Didinium nasutum and three Paramecium sp, which differ in pulsed-field gel electrophoresis (PFGE) karyotype, and compared with the model of topologically associated domains (TADs) of higher eukaryotic nuclei. PFGE showed that the sizes of macronuclear DNAs ranged from 50 to 1700 kb, while the majority of the molecules were less than 500 kb in length. A comparative quantitative analysis of the PFGE and electron microscopic data showed that each chromatin body contained one minichromosome in P. multimicronucleatum in the logarithmic growth phase, while bodies in the D. nasutum macronucleus contained two or more DNA molecules each. Chromatin bodies aggregated during starvation, when activity of the macronuclei decreased, leading to an increase of chromatin body size or the formation of 200- to 300-nm fibrils of several chromatin bodies. A model was proposed to explain the formation of such structures. In terms of topological characteristics, macronuclear chromatin bodies with subchromosomal DNA molecules were found to correspond to higher eukaryotic TADs.

Macronuclear Plasticity in Two South American Populations of Spirostomum (Ciliophora, Heterotrichea) Warns About its Use for Species Classification: Revision and New Insights.

Spirostomum is a widely distributed heterotrichean genus composed of well-known species with described ecology and phylogenetic affinities. The morphological classification of Spirostomum species is mostly based on the body size/shape, number of cortical granule rows and macronuclear characteristics. These features along with molecular phylogenies based on ribosomal genes divide the genus into two phylogroups, one including species with a compact macronucleus, and another including species with a moniliform macronucleus. Here, we present our observations on atypical Spirostomum specimens with unusually two distinct macronuclei and shortened adoral zone of membranelles. These atypical forms appeared in the cultures of S. minus and S. yagiui, sampled at different sites in South America (Chile and Brazil) and associated with unrelated substrate types. Morphological observations of living and stained cells, 18S rRNA gene analyses, and a thorough investigation of the literature suggest that the atypical phenotype may be a result of uncommon pathways during the conjugative process. Thus, we demonstrate that studies of ciliate natural populations and their morphological variations, especially from undersampled biogeographical regions, can reveal the boundaries of widely used morphological characters for Spirostomum taxonomy and species identification.

Characterization and use of Tetrahymena thermophila artificial chromosome 2 (TtAC2) constructed by biomimetic of macronuclear rDNA minichromosome.

Efficient expression vectors for unicellular ciliate eukaryotic Tetrahymena thermophila are still needed in recombinant biology and biotechnology applications. Previously, the construction of the T. thermophila Macronuclear Artificial Chromosome 1 (TtAC1) vector revealed additional needs for structural improvements such as better in vivo stability and maintenance as a recombinant protein expression platform. In this study, we designed an efficiently maintained artificial chromosome by biomimetic of the native macronuclear rDNA minichromosome. TtAC2 was constructed by sequential cloning of subtelomeric 3'NTS region (1.8 kb), an antibiotic resistance gene cassette (2 kb neo4), a gene expression cassette (2 kb TtsfGFP), rDNA coding regions plus a dominant C3 origin sequence (10.3 kb), and telomeres (2.4 kb) in a pUC19 backbone plasmid (2.6 kb). The 21 kb TtAC2 was characterized using fluorescence microscopy, qPCR, western blot and Southern blot after its transformation to vegetative T. thermophila CU428.2 strain, which has a recessive B origin allele. All experimental data show that circular or linear forms of novel TtAC2 were maintained as free replicons in T. thermophila macronucleus with or without antibiotic treatment. Notably, TtAC2 carrying strains expressed a TtsfGFP marker protein, demonstrating the efficacy and functionality of the protein expression platform. We show that TtAC2 is functionally maintained for more than two months, and can be efficiently used in recombinant DNA, and protein production applications.

The Compact Macronuclear Genome of the Ciliate Halteria grandinella: A Transcriptome-Like Genome with 23,000 Nanochromosomes.

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions. This challenges our usual understanding of chromosomal structure and suggests the possibility of novel mechvanisms in transcriptional regulation. Comprehensive analysis of multiple data sets reveals that Halteria transcription dynamics are influenced by: (i) nonuniform nanochromosome copy numbers correlated with gene-expression level; (ii) dynamic alterations at both the DNA and RNA levels, including alternative internal eliminated sequence (IES) deletions during macronucleus formation and large-scale alternative splicing in transcript maturation; and (iii) extremely short 5' and 3' untranslated regions (UTRs) and universal TATA box-like motifs in the compact 5' subtelomeric regions of most chromosomes. This study broadens the view of ciliate biology and the evolution of unicellular eukaryotes, and identifies Halteria as one of the most compact known eukaryotic genomes, indicating that complex cell structure does not require complex gene architecture.

The Function and Evolution of Motile DNA Replication Systems in Ciliates.

DNA replication is a ubiquitous and conserved cellular process. However, regulation of DNA replication is only understood in a small fraction of organisms that poorly represent the diversity of genetic systems in nature. Here we used computational and experimental approaches to examine the function and evolution of one such system, the replication band (RB) in spirotrich ciliates, which is a localized, motile hub that traverses the macronucleus while replicating DNA. We show that the RB can take unique forms in different species, from polar bands to a "replication envelope," where replication initiates at the nuclear periphery before advancing inward. Furthermore, we identify genes involved in cellular transport, including calcium transporters and cytoskeletal regulators, that are associated with the RB and may be involved in its function and translocation. These findings highlight the evolution and diversity of DNA replication systems and provide insights into the regulation of nuclear organization and processes.

Identification and Characterization of Base-Substitution Mutations in the Macronuclear Genome of the Ciliate Tetrahymena thermophila.

Polyploidy can provide adaptive advantages and drive evolution. Amitotic division of the polyploid macronucleus (MAC) in ciliates acts as a nonsexual genetic mechanism to enhance adaptation to stress conditions and thus provides a unique model to investigate the evolutionary role of polyploidy. Mutation is the primary source of the variation responsible for evolution and adaptation; however, to date, de novo mutations that occur in ciliate MAC genomes during these processes have not been characterized and their biological impacts are undefined. Here, we carried out long-term evolution experiments to directly explore de novo MAC mutations and their molecular features in the model ciliate, Tetrahymena thermophila. A simple but effective method was established to detect base-substitution mutations in evolving populations whereas filtering out most of the false positive base-substitutions caused by repetitive sequences and the programmed genome rearrangements. The detected mutations were rigorously validated using the MassARRAY system. Validated mutations showed a strong G/C→A/T bias, consistent with observations in other species. Moreover, a progressive increase in growth rate of the evolving populations suggested that some of these mutations might be responsible for cell fitness. The established mutation identification and validation methods will be an invaluable resource to make ciliates an important model system to study the role of polyploidy in evolution.

Genome plasticity in Paramecium bursaria revealed by population genomics.

BACKGROUND: Ciliates are an ancient and diverse eukaryotic group found in various environments. A unique feature of ciliates is their nuclear dimorphism, by which two types of nuclei, the diploid germline micronucleus (MIC) and polyploidy somatic macronucleus (MAC), are present in the same cytoplasm and serve different functions. During each sexual cycle, ciliates develop a new macronucleus in which newly fused genomes are extensively rearranged to generate functional minichromosomes. Interestingly, each ciliate species seems to have its way of processing genomes, providing a diversity of resources for studying genome plasticity and its regulation. Here, we sequenced and analyzed the macronuclear genome of different strains of Paramecium bursaria, a highly divergent species of the genus Paramecium which can stably establish endosymbioses with green algae. RESULTS: We assembled a high-quality macronuclear genome of P. bursaria and further refined genome annotation by comparing population genomic data. We identified several species-specific expansions in protein families and gene lineages that are potentially associated with endosymbiosis. Moreover, we observed an intensive chromosome breakage pattern that occurred during or shortly after sexual reproduction and contributed to highly variable gene dosage throughout the genome. However, patterns of copy number variation were highly correlated among genetically divergent strains, suggesting that copy number is adjusted by some regulatory mechanisms or natural selection. Further analysis showed that genes with low copy number variation among populations tended to function in basic cellular pathways, whereas highly variable genes were enriched in environmental response pathways. CONCLUSIONS: We report programmed DNA rearrangements in the P. bursaria macronuclear genome that allow cells to adjust gene copy number globally according to individual gene functions. Our results suggest that large-scale gene copy number variation may represent an ancient mechanism for cells to adapt to different environments.

The completed macronuclear genome of a model ciliate Tetrahymena thermophila and its application in genome scrambling and copy number analyses.

The ciliate Tetrahymena thermophila has been a powerful model system for molecular and cellular biology. However, some investigations have been limited due to the incomplete closure and sequencing of the macronuclear genome assembly, which for many years has been stalled at 1,158 scaffolds, with large sections of unknown sequences (available in Tetrahymena Genome Database, TGD, http://ciliate.org/ ). Here we completed the first chromosome-level Tetrahymena macronuclear genome assembly, with approximately 300× long Single Molecule, Real-Time reads of the wild-type SB210 cells-the reference strain for the initial macronuclear genome sequencing project. All 181 chromosomes were capped with two telomeres and gaps were entirely closed. The completed genome shows significant improvements over the current assembly (TGD 2014) in both chromosome structure and sequence integrity. The majority of previously identified gene models shown in TGD were retained, with the addition of 36 new genes and 883 genes with modified gene models. The new genome and annotation were incorporated into TGD. This new genome allows for pursuit in some underexplored areas that were far more challenging previously; two of them, genome scrambling and chromosomal copy number, were investigated in this study. We expect that the completed macronuclear genome will facilitate many studies in Tetrahymena biology, as well as multiple lines of research in other eukaryotes.

Functional diversification of Paramecium Ku80 paralogs safeguards genome integrity during precise programmed DNA elimination.

Gene duplication and diversification drive the emergence of novel functions during evolution. Because of whole genome duplications, ciliates from the Paramecium aurelia group constitute a remarkable system to study the evolutionary fate of duplicated genes. Paramecium species harbor two types of nuclei: a germline micronucleus (MIC) and a somatic macronucleus (MAC) that forms from the MIC at each sexual cycle. During MAC development, ~45,000 germline Internal Eliminated Sequences (IES) are excised precisely from the genome through a 'cut-and-close' mechanism. Here, we have studied the P. tetraurelia paralogs of KU80, which encode a key DNA double-strand break repair factor involved in non-homologous end joining. The three KU80 genes have different transcription patterns, KU80a and KU80b being constitutively expressed, while KU80c is specifically induced during MAC development. Immunofluorescence microscopy and high-throughput DNA sequencing revealed that Ku80c stably anchors the PiggyMac (Pgm) endonuclease in the developing MAC and is essential for IES excision genome-wide, providing a molecular explanation for the previously reported Ku-dependent licensing of DNA cleavage at IES ends. Expressing Ku80a under KU80c transcription signals failed to complement a depletion of endogenous Ku80c, indicating that the two paralogous proteins have distinct properties. Domain-swap experiments identified the α/ß domain of Ku80c as the major determinant for its specialized function, while its C-terminal part is required for excision of only a small subset of IESs located in IES-dense regions. We conclude that Ku80c has acquired the ability to license Pgm-dependent DNA cleavage, securing precise DNA elimination during programmed rearrangements. The present study thus provides novel evidence for functional diversification of genes issued from a whole-genome duplication.

Rearrangement of macronucleus chromosomes correspond to TAD-like structures of micronucleus chromosomes in Tetrahymena thermophila.

The somatic macronucleus (MAC) and germline micronucleus (MIC) of Tetrahymena thermophila differ in chromosome numbers, sizes, functions, transcriptional activities, and cohesin complex location. However, the higher-order chromatin organization in T. thermophila is still largely unknown. Here, we explored the higher-order chromatin organization in the two distinct nuclei of T. thermophila using the Hi-C and HiChIP methods. We found that the meiotic crescent MIC has a specific chromosome interaction pattern, with all the telomeres or centromeres on the five MIC chromosomes clustering together, respectively, which is also helpful to identify the midpoints of centromeres in the MIC. We revealed that the MAC chromosomes lack A/B compartments, topologically associating domains (TADs), and chromatin loops. The MIC chromosomes have TAD-like structures but not A/B compartments and chromatin loops. The boundaries of the TAD-like structures in the MIC are highly consistent with the chromatin breakage sequence (CBS) sites, suggesting that each TAD-like structure of the MIC chromosomes develops into one MAC chromosome during MAC development, which provides a mechanism of the formation of MAC chromosomes during conjugation. Overall, we demonstrated the distinct higher-order chromatin organization in the two nuclei of the T. thermophila and suggest that the higher-order chromatin structures may play important roles during the development of the MAC chromosomes.

Deciphering phylogenetic relationships and delimiting species boundaries using a Bayesian coalescent approach in protists: A case study of the ciliate genus Spirostomum (Ciliophora, Heterotrichea).

The ciliate genus Spirostomum comprises eight morphospecies, inhabiting diverse aquatic environments worldwide, where they can be used as water quality indicators. Although Spirostomum species are relatively easily identified using morphological methods, the previous nuclear rDNA-based phylogenies indicated several conflicts in morphospecies delineation. Moreover, the single locus phylogenies and previous analytical approaches could not unambiguously resolve phylogenetic relationships among Spirostomum morphospecies. Here, we attempt to investigate species boundaries and evolutionary history of Spirostomum taxa, using 166 new sequences from multiple populations employing one mitochondrial locus (CO1 gene) and two nuclear loci (rRNA operon and alpha-tubulin gene). In accordance with previous studies, relationships among the eight Spirostomum morphospecies were poorly supported statistically in individual gene trees. To overcome this problem, we utilised for the first time in ciliates the Bayesian coalescent approach, which accounts for ancestral polymorphisms, incomplete lineage sorting, and recombination. This strategy enabled us to robustly resolve deep relationships between Spirostomum species and to support the hypothesis that taxa with compact macronucleus and taxa with moniliform macronucleus each form a distinct lineage. Bayesian coalescent-based delimitation analyses strongly statistically supported the traditional morphospecies concept but also indicated that there are two S. minus-like cryptic species and S. teres is non-monophyletic. Spirostomum teres was very likely defined by a set of ancestral features of lineages that also gave rise to S. yagiui and S. dharwarensis. However, molecular data from type populations of the morphospecies S. minus and S. teres are required to unambiguously resolve the taxonomic problems.

Tetrahymena Comparative Genomics Database (TCGD): a community resource for Tetrahymena.

Ciliates are a large and diverse group of unicellular organisms characterized by having the following two distinct type of nuclei within a single cell: micronucleus (MIC) and macronucleus (MAC). Although the genomes of several ciliates in different groups have been sequenced, comparative genomics data for multiple species within a ciliate genus are not yet available. Here we collected the genome information and comparative genomics analysis results for 10 species in the Tetrahymena genus, including the previously sequenced model organism Tetrahymena thermophila and 9 newly sequenced species, and constructed a genus-level comparative analysis platform, the Tetrahymena Comparative Genomics Database (TCGD). Genome sequences, transcriptomic data, gene models, functional annotation, ortholog groups and synteny maps were built into this database and a user-friendly interface was developed for searching, visualizing and analyzing these data. In summary, the TCGD (http://ciliate.ihb.ac.cn) will be an important and useful resource for the ciliate research community.